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ABSTRACT 

Previous submillimetre (submm) observations detected 0.7 M Q of cool dust emission around 
the Luminous Blue Variable (LBV) star i] Carinae. These observations were hindered by the 
low declination of r\ Carinae and contamination from free-free emission orginating from the 
stellar wind. Here, we present deep submm observations with LABoCa at 870 /im, taken 
shortly after a maximum in the 5.5-yr radio cycle. We find a significant difference in the 
submm flux measured here compared with the previous measurement: the first indication of 
variability at submm wavelengths. A comparison of the submm structures with ionised emis- 
sion features suggests the 870 /im is dominated by emission from the ionised wind and not 
thermal emission from dust. We estimate 0.4±0.1 M Q of dust surrounding 77 Carinae. The spa- 
tial distribution of the submm emission limits the mass loss to within the last thousand years, 
and is associated with mass ejected during the great eruptions and the pre-outburst LBV wind 
phase; we estimate that T] Carinae has ejected > 40 M of gas within this timescale. 
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1 INTRODUCTION 

77 Carinae (77 Car), one of the most luminous infra-red (IR) objects 
in our Galaxy (>10 6 Lq), is well known for dramatic outbursts over 
the last thousands of years, in which vast amounts of material is 
ejected outwards from the star. The great eruption phase during the 
last 200 years are responsible for creating the dusty bipolar nebula 
known as the Homunculus, as seen in the Hubble Space Telescope 
(HST) images. The recent discovery of a fast blast wave from the 
1843 eruption (Smith 2008) indicates that the mass loss in 77 Car 
is not simply due to a stellar wind, but could be more akin to a 
low-energy supernova remnant. Extensive IR imaging has revealed 
important clues about the mass loss history; the Infrared Space Ob- 
servatory (ISO) discovered a massive cool dust component (Morris 
et al. 1999; Smith et al. 2003) with a total of ~O.1M of dust 
proposed to be present in the nebula, suggesting a mass loss rate of 
0.5 Mq yr _1 during the last 200 years. Submillimetre (submm) ob- 
servations (Gomez et al. 2006, hereafter G06) of 77 Car at 450 and 
850 77m with the Submillimetre Common User Bolometer Array 
(SCUBA) detected the presence of a massive component of dust, 
with ~O.7M0 needed to reproduce the IR-submm Spectral En- 
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ergy Distribution (SED). Their work suggested that up to four times 
more mass had been ejected from 77 Car during its recent, violent 
history than previously proposed. The SCUBA data also indicated 
that the cool dust component was extended along the mid-plane, 
far beyond the inner optical and IR region; this was interpreted as 
mass which had been ejected on a much longer timescale than the 
well-known, smaller features such as the torus and Homunculus. 

The previous analysis of the SCUBA data was hindered by 
two unknowns: whether the submm emission originated from free- 
free and not thermal emission from dust and whether the extended 
structure was reliable. Firstly, the submm fluxes are likely to be 
contaminated by strong free-free emission from the extended stel- 
lar wind which is thought to vary with frequency as j, - 6-1 - 3 ( e g p 
Cox et al. 1995; Brooks et al. 2005, hereafter B05). This is fur- 
ther compounded since the radio and millimetre (mm) fluxes are 
highly variable, with 7/ Car at 3 cm changing appearance from a 
point source in 1992 to an extended region with five times more 
flux in 1996 (Duncan & White 2003). The X-ray, radio and mm 
variability of 77 Car has been well documented (Pittard et al. 1998; 
Abraham et al. 2005; Damineli et al. 2008), with the star under- 
going periodic variability over a 5.5-yr period. This variability has 
been associated with shocks from an extended disc surrounding the 
77 Car binary system (e.g. Duncan & White 2003) and/or due to 
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episodic mass loss which, in turn, decreases the number of ionising 
photons. Since there were no available mm observations to match 
the epoch of the original 1998 SCUBA data (taken during the min- 
imum phase of the radio cycle), G06 were not able to estimate the 
contribution to the submm due to free-free emission. Second, the 
SCUBA data were taken at very large airmass (A ~ 5.5) due to the 
declination of r\ Car. This meant that the beam shape at 450 pm was 
distorted due to the deformation of the dish at such low elevations 
and the morphology of the source was difficult to ascertain from 
this dataset. Accurately determining the distribution and mass of 
dust in the stellar wind has enormous consequences not only on de- 
termining the mass loss history of r\ Car, but also understanding the 
progenitors of 'exotic' core-collapse supernovae (SNe; Pastorello 
et al. 2007; Smith et al. 2009). 

In this Letter, we present deep submm LABoCa observations 
of r\ Car, taken shortly after the maximum phase in the cycle. In ij2] 
we discuss the data reduction and compare the 870 pm emission to 
well-known features seen in X-ray, optical, Ha, 8 /im and 1.2 mm 
in S|3] In Sj4]we estimate the dust and gas mass for the stellar wind 
and our conclusions are presented in §J3] 



2 OBSERVATIONS AND DATA REDUCTION 

The 870 pm data were taken with the LABoCa camera (Siringo 
et al. 2009), a 295-pixel bolometer array, located on the Atacama 
Pathfinder Experiment telescope (Giisten et al. 2006) on Chanjan- 
tor in Chile. r\ Car was observed during Science Verification on 
23rd July 2007. This epoch corresponds to 1.4 yrs after maximum 
brightness in the radio and X-ray cycles (2006.25). The observa- 
tions were carried out in raster spiral mode with twelve scans, each 
providing a fully sampled map over 30' x 30'. The total on-source 
integration time was 1 .7 hours. Two independant mesurements of 
the optical depth, r were obtained. The first method used the pre- 
cipitable water vapour (PWV) levels measured every minute along 
the line of sight, then scaled using the relevant atmopsheric trans- 
mission model. The PWV ranged from 0.7-0.9 mm. The second 
method calculated r from skydip measurements, where a model of 
the dependence of the effective sky temperature on elevation were 
fitted to determine the zenith opacity. The two skydips taken before 
and after the on-source scans were both well fitted by the model, 
with r between 0.1-0.2. These values were 25% lower than those 
estimated from the PWV measurements. A linear combination of 
the two methods, and a final comparison with the calibrator models 
(e.g. Siringo et al. 2009), was used to determine the values used 
in the data reduction. Following Dunne & Eales (2001), the frac- 
tional error in the flux from the error due to the range of opacities 
calculated with these two methods is 6 per cent. 

The data were reduced using the BoA (BOlometer array Anal- 
ysis software) package. The focus was checked on observations of 
Venus and Jupiter and was stable within ±0.2 mm. Pointing obser- 
vations were within 1" in azimuth and elevation during the earlier 
scans but crept up to within 3" towards the end. Bad and noisy 
pixels were flagged, with 276 bolometers used in each scan. The 
data were despiked and correlated noise was removed. The reduc- 
tion is optimized for the recovery of strong sources. The scans were 
coadded (weighted by rms -2 ) and the data were gridded onto 6" 
pixels to create the final map. After a first iteration of the reduc- 
tion, the source map was used to flag bright sources and the data 
were reduced again. This was efficient at removing negative arti- 
facts which appear around the bright sources in the first iteration 
and led to a more stable background noise level in the central re- 



gion (~40mJy beam -1 ). The 870 pm map towards r\ Car is shown 
in Fig. Q] together with 8 pm Midcourse Space experiment (MSX) 
archival data. r\ Car is the strongest source in the LABoCa map, 
and peaks at signal-to-noise >100. Well known emission regions, 
molecular clouds and star clusters are identified in the image. 

Secondary calibrators were used to calibrate the map and the 
typical correction factor when comparing the measured flux on the 
calibrator with its expected flux in an aperture was ~5 per cent. The 
total uncertainty in the calibration, including uncertainties in the 
calibrator model is therefore ~ 12 per cent. The 870 pm flux of r\ 
Car measured in an aperture with radius 30" is 42 ± 5 Jy, compared 
to the flux estimated by SCUBA in 1998 of 13 ± 1 Jy (G06). This 
is an increase by a factor of three over the 9 yr period and cannot 
be accounted for by calibration errors. 



3 RESULTS 

A zoomed in image of r\ Car is presented in Fig. Q] and compared 
with Ha (Smith, Bally & Morse. 2003) and 1.2 mm (B05) emis- 
sion. Faint structures not previously seen in FIR-submm maps are 
detected with LABoCa and are identified with well known sources 
(labeled in Fig. Q] following the notation of Smith et al. 2003 & 
B05). The submm emission can be separated into four components: 
the bright source at the location of r\ Car; filaments and arcs which 
surround molecular globules (Rathbome et al. 2002); peaks and fil- 
aments associated with molecular clouds at the edges of the map; 
and a number of compact sources (seen at mm and submm wave- 
lengths, but not always in Ha and 8 pm). The edges of the two 
large molecular clouds, part of a giant cloud complex within the 
Carina arm, are visible in the LABoCa map to the North and South 
of r\ Car. The 870 pm flux in these regions correlates with peaks 
seen in molecular emission (Cox et al. 1995; Yonekura et al. 2005) 
suggesting dust emission from cold clouds. 

r/ Car is the brightest source at 870 pm and 1.2 mm, but not at 
IR-FIR. At 12 and 25 pm, the peak is found at the centre of Tr 14 
(a massive star cluster) whereas at 60 and 100 pm, the peak shifts 
towards the west (and the Car I feature). This shift traces a tempera- 
ture gradient in the nebula (Cox 1995); hotter dust is located closer 
to Tr 14, where massive stars are externally heating the clouds and 
cooler material is forced further away. The shift in peak emission 
from IR-FIR is also seen in our submm data. 

The faint structure just north of 77 Car is associated with the 
Car II emission feature. The 870 pm emission closely follows the 
Ha, 8 pm and 1.2 mm emission (Fig.[2j» and at longer wavelengths, 
the 3 cm radio emission (B05). These arcs and filaments surround 
dark clouds detected in CO emission (Rathbome et al. 2002, Cox & 
Bronfman 1995). The close correlation between ionised gas, emis- 
sion at both mm and radio wavelengths and the LABoCa emis- 
sion suggest that the 870 pm originates from ionised material. In- 
deed, the integrated flux of Car II at 870 pm here agrees with 
the extrapolated radio continuum spectrum in B05 (their Fig 5) 
which are well fitted with an Hll emission model. The ratio of 
the 870 pm/1.2 mm maps confirm that the spectral index in Car II 
varies as F v tx v~ 2 , consistent with an ionization front or Hll 
region, thus the LABoCa emission here orginates from a free-free 
source and not from warm dust. Our lack of submm detection at the 
location of the famous Keyhole Nebula (labeled 1-3, Fig. [TJ rules 
out any dust component with temperatures higher than 13 K. 
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Figure 1. The environment surrounding rj Car. Top Panel: Left: LABoCa 870 fim emission. Right: MSX 8 fim emission, shown here in negative log greyscale, 
with LABoCa 870 fim contours (starting at 0.04Jybeam _1 ). Region shown covers 19' X 19' area with LABoCa beam size shown in lower right hand side. 
Known compact sources are labeled 4-8 (B05). Bottom Panel: Zoomed in region (9' X 8.4') with LABoCa 870 fim contours overlaid, ranging from 0.04- 
0.25 Jy beam" 1 . Left. LABoCa emission. Middle: 1.2 mm emission (SIMBA map kindly provided by K Brooks B05). Right: Ha emission shown here in 
negative greyscale (CTIO map kindly provided by N Smith, Smith, Bally & Morse 2003). The Keyhole Nebula is labeled 1, 2 and 3. 



3.1 Submillimetre Emission from r\ Carinae 

The central r\ car source as seen with LABoCa (Figs. [T] & O is 
well fitted by a two-dimensional Gaussian with FWHM 22". We 
find that r\ car is unresolved with the LABoCa beam (19") and 
we see the first error beam pattern below 1 per cent peak intensity 
at 0.3 Jybeam -1 . Deconvolving from the beam and pixel size in- 
dicates a central source with FWHM ~10". We do not detect the 
structure along the midplane seen in the original 850 fim data from 
G06, confirming the suggestion made there that a combination of 
high airmass and chop throw had smeared the signal across the mid- 
plane. Based on the LABoCa observations, the true extent of the 
submm source is likely to be within the Homunculus (18") and in- 
terior to the region encompassed by the O III 'cocoon' of material 
(60") detected by Smith, Morse & Bally (2005). Recently, Smith 
(2008) detected a fast blast wave, traveling at up to 6000 km s , 
from the great eruption in 1843. This blast wave extends out to the 
same radius as the cocoon, and is interior to or coincident with, 
the X-ray emission. Smith proposed that the X-rays arise from this 
fast blast wave running into a previous eruption which gave rise to 
a nitrogen rich shell. Such a scenario places severe constraints on 



the timescales for emission seen interior to the blast wave material 
and suggests that the source responsible for the emission seen by 
LABoCa was ejected within the last 200-1000 years. 

The 1 .2 mm emission from r\ Car is closely correlated with the 
870 /im, and neither correlate well with the IR-FIR central source. 
If we assume that the former originates from thermal emission 
from dust grains (i.e. F v oc v' i+l3 ), the ratio of 870/im/1.2mm 
would require a dust emissivity index, 0, of -2. If we assume emis- 
sion at both wavelengths originates from a source which varies as 
F v oc v~ a , we obtain a — 0.5-2 which agrees well with emission 
from an ionised stellar wind source. The similar variability seen at 
both 870 /im and 1.2 mm suggests that less than 20 per cent of the 
870 /im flux is contributed by dust emission. 



4 DISCUSSION 

4.1 The Spectral Energy Distribution of rj Car 

The spectral energy distribution (SED) of r\ Car is shown in Fig. [3] 
Given the variability of observed fluxes over different epochs dur- 
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Figure 2. LABoCa 870 fim emission of r\ Car with archival X-ray (CHAN- 
DRA) and optical (HST) contours shown in black and white respectively. 
The LABoCa beam (19") is shown in the lower right hand side. Units are 
Jy beam - 1 . The extent of the Homunculus seen as the double-lobe structure 
in the contours is ~ 18". 
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Figure 3. The SED of r\ Car from 12 fim through to the radio. The radio, 
mm and 850/870 fim fluxes are shown as vertical bars (following B05) to 
indicate their variability range during the 5.5-yr cycle. The solid black line 
is the sum of the two-temperature components fitted to the IR (T\ = 174 K 
and T w = 82 K) and 450 fim fluxes (dot-dashed and dashed respectively). 
The dot-dot-dot-dashed line shows the SED expected from the optically 
thick free-free wind, F v oc u 6 and the optically thin wind with oc u l s . 
The dashed line represents emission from the optically thin Little Ho- 
munculus feature tx v~ . The red line shows the total SED from all the 
components. The far-IR luminosity from 12-1000 fim is 1.6 X 10 6 Lq. 



ing rj Car's radio cycle, we represent the variability as vertical bars 
following B05, highlighting the flux changes over the cycle rather 
than the fluxes measured at different epochs. 

There are three major components to the SED: thermal emis- 
sion from dust grains (following a modified blackbody with power- 
law index /3), free-free emission from the ionised stellar wind and 
free-free emission from the optically thin Little Homunculus fea- 
ture, seen close to the central star (Duncan & White 2003). The 
stellar wind material consists of two separate components: the op- 
tically thin varying as u 13 (e.g. B05) and the optically thick vary- 
ing as a classical wind with v ' 6 (Wright & Barlow 1975; Lamers 
& Cassinelli 1999). The optically thin component from the Ho- 
munculus is produced by the ionisation of rj Car's stellar wind by 
the UV-radiation field of the hot binary companion. This produces 
free-free emission which varies as u~ 0A and relates to the Little 
Homunculus and torus (Teodoro et al. 2008). 

For the dust component, we used two-modified blackbodies: 

Su = N h v^B{u,T h )+N w v p B{u,T w ) (1) 

where N is the normalisation term, (3 is the dust emissivity 
index and B is the Planck function at frequency v and tempera- 
ture T. The SCUBA 850 and LABoCa 870 fim fluxes are not in- 
cluded in the fit to the thermal emission from dust. The 450 fim 
was corrected for free-free emission expected at this wavelength 
using the predicted power-law variation u 6 (as shown in Fig. [3) 
where ~35 per cent of the 450 1998 flux could originate from 
free-free emission. The revised 450 fim flux due to dust in the stel- 
lar wind is therefore ~30 Jy. This agrees with the results from the 
Balloon-borne Large Aperture Submillimeter Telescope who ob- 
served rf Car from 250-500 fim during the maximum in 2006 (Har- 
grave priv. comm., Hargrave et al. in prep). Hargrave et al. find 
no evidence of an increased flux at 500 fim when compared to the 
SCUBA 450 fim observed during the minimum. This supports the 
idea that the 450 fim flux measured with SCUBA is still dominated 



by thermal emission from dust grains and does not undergo signif- 
icant variability. Note that the 450 flux is key to the dust mass 
estimation since this extra constraint forces the fit to have a lower 
dust temperature and a steeper emissivity index than the previous 
published models (Morris et al. 1999; Smith et al. 2003; B05). We 
rule out a cold dust component since we cannot obtain a modi- 
fied blackbody fit to the submm-mm fluxes. An adequate fit can be 
obtained using only the mm fluxes but this would require a dust 
component with T~4 K; it is difficult to understand how such cold 
dust could survive the close environment of rf Car. 
The dust mass is estimated using Eq.[2] 

M50 \B(U450,lh) B(U450,J-w)J 

where k„ is the dust-mass-absorption coefficient 
(0.27kg _1 m 2 for interstellar dust at 450 fim, James et al. 
2002) and d is the distance to the star. The best-fit (\ ) to the SED 
requires a total dust mass of 0.4 ± 0.1 Mq and /3~1.9 (indicative 
of normal interstellar dust grains), T h =174^ 6 K and T w =80± 1° K. 
The errors are estimated from the bootstrap technique where 1000 
fits to the fluxes were made; the 68% confidence intervals are 
quoted here. Previous authors required /3 ~ 1 to fit the mid-IR 
SED but higher values are necessary to fit the long-wavelength 
submm fluxes at 450 and 870 fim. 

The "transition" between the classical u ' 6 and i/ 1,3 power- 
law regimes from the stellar wind as plotted in Fig. [3] was noted in 
Cox et al. (1995) and B05 and is required to explain the observed 
SED at mm wavelengths. They proposed that this transition occurs 
due to changes in the ionisation properties of the wind but the lo- 
cation of this transition is not well constrained and consequently 
affects the amount of free-free emission expected at submm wave- 
lengths. For example, if the two power-law regimes turned-over at 
870 instead of 1.2 mm as plotted here, this would suggest that 
100 per cent of the 870 fim flux and 65 per cent of the 450 could 
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originate from a free-free source. Changes to the wavelength of the 
transition will change the derived properties which fit the thermal 
dust SED (T, P) and the dust mass, but these changes are well 
within the errors quoted above. The exact location of this transi- 
tion, though important to understand for physical reasons, does not 
affect the dust mass estimated in the current work. 

With gas-to-dust ratio ~100 (Smith & Ferlund 2007), we es- 
timate that >40 M Q of gas has been ejected. This is approximately 
half of the value quoted in G06 but three times more than the max- 
imum values quoted from previous IR data (Morris et al. 1999; 
Smith et al. 2003). Our estimate is also consistent with the up- 
per end of the mass range (15-35 M ) found by Smith & Ferlund 
(2007) using an independent, theoretical model, which determined 
the gas mass within the Homunculus using the recent detection of 
H2 to determine the density of the gas. Note that the dust mass 
quoted here is different to G06, since the submm flux attribruted to 
dust emission has decreased in this work, yet the dust temperature 
required to fit the steeper slope between 175 (im and 1000 /im is 
colder here, thereby increasing the mass. 

4.2 Mass Loss traced by the Submillimetre Emission 

The LABoCa submm emission is interior to the X-ray shell (tracing 
blast wave material), the O III veil (tracing the recent N-rich mass 
loss interacting with older O-rich material (Smith et al. 2005) and 
the fast blast wave from recent ejecta interacting with the older N- 
rich material (Smith 2008). Since we have little information about 
the true structure of the submm source, the LABoCa emission is 
most likely a complex combination of the mass lost during the 
great ejection phases over the last few hundred years, and the pre- 
outburst stellar wind material ejected over the last few thousand 
years. From our constraints on the spatial location of the dust, we 
therefore estimate an average mass loss rate of >10~ 2 M Q yr _1 , 
at least one order of magnitude higher than suggested in Hillier et 
al. (2001) and Smith et al. (2005). We are currently unable to re- 
solve the features within the Homunculus and pin down the exact 
distribution or contributions of the free-free and dust components, 
hence our mass loss has a large associated uncertainty. We do not 
know the exact timescale over which the mass has been lost nor 
the composition and radiation properties of the dust grains. The 
longer timescale put on the submm mass loss in this work com- 
pared to G06, eases the strain on the available mass reservoir of r\ 
Car (Smith 2009b), since the 40 M Q of gas traced by the submm 
is likely a combination of eruptions and the stellar wind. However, 
we note that such a huge mass loss does provide support for the 
idea of dredge-up SN-like events and also hints that pre-SN mass 
loss could be significant contributors to the interstellar dust budget 
in the early Universe (Morgan & Edmunds 2003; Dwek, Galliano 
& Jones 2007). 



5 CONCLUSIONS 

We found that the submm extended emission surrounding r\ Car, 
originally identified in G06, is unresolved by the LABoCa beam 
at 870 /im, and is highly variable, following the well-documented 
cycle seen at mm and radio wavelengths. Accounting for the flux 
difference from G06 and associating the 800-1000 /im with free- 
free emission rather than thermal emission from grains, we revise 
the dust mass in the stellar wind to 0.4 ±0.1 M Q . This translates 
to a mass loss greater than 40 M© over the past 1000 years. Future 
observations with Herschel and ALMA are needed to constrain the 



submm and mm variability and multi-epoch data will allow us to 
separate out the dust and free-free components throughout the cy- 
cle. Higher angular resolution data are needed to separate out the 
multiple phases of mass loss, this will provide crucial information 
about the evolutionary path to massive-star SN explosions. 
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